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Abstract 


This paper introduces the measure of approximate-degree and the concept of 
approximate-degree function between numerical values, thus developing a new 
interpolation method —— approximation-degree-based interpolation, i.e., AD 
interpolation. One-dimensional AD interpolation is done directly by using correlative 
interpolation formulas; n(n>1)-dimensional AD interpolation is firstly separated into n 
parallel one-dimensional AD interpolation computations to do respectively, and then 
got results are synthesized by Sum-Times-Difference formula into a value as the result 
value of the n-dimensional interpolation. If the parallel processing is used, the 
efficiency of n-dimensional AD interpolation is almost the same as that of the 
one-dimensional AD interpolation. Thus it starts a feasible and convenient approach 
and provides an effective method for high-dimensional interpolations. Furthermore, if 
AD interpolation is introduced into machine learning, a new instance-based learning 


method is expected to be realized. 


Keywords Approximation-Degree; High-Dimensional Interpolation; First Separating 
then Synthesizing; Parallel Processing; Sum-Times-Difference Formula; 


Instance-Based Learning 


1 Introduction 


With the rapid development of artificial intelligence and machine learning, some of the 
traditional numerical computation methods such as interpolation, data fitting and regression 
analysis have aroused people's interest and become active fields again. Although there has 
been sufficient research on interpolation, the results of high-dimensional interpolation are still 
relatively few, and some of the existing high-dimensional interpolation methods"! have a large 
load of computation and are not suitable for local interpolation. High-dimensional 
interpolation has important applications in artificial intelligence and machine learning. For 


BA often involves high-dimensional interpolation. Therefore, 


example, instance-based learning 
the high-dimensional interpolation, especially that suitable for local interpolation, is still an 
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important subject worthy of serious study. 

Inspired by the approximate evaluation method of flexible linguistic functions ° in 
reference [1], this paper intends to introduce the measure of approximate-degree between 
numerical values to study the approximate evaluation of numerical functions, and then 
explores the approach and method of interpolation based on approximate-degree, especially, 


the method of high-dimensional interpolation based on approximate-degree. 


2 Approximation Axiom, Approximation-Degree, and 


Approximation-Degree Transmission 


Approximation Axiom Let R be real number field, R” be n-dimensional real vector 
space, U=[a1, bı]x[az, b2]x...X[dn, bn] ER” (n=1),V=[c, d]cR, and y=f(x) (x=(1, X2,...5 Xn)) 
be a continuous and non-chaotic function relation from U to V. If yo= f(xo) (xoEU) is known, 
then when xeU is approximate to xo, y= f(x) €V is also approximate to yo. 

Definition 2-1 Let R be real number field, xoe[a, b]cR, and [a@, folc[a, b] be a 
neighborhood of xo, which is called the approximation domain of xo. For Vxe[a, b], say x is 
approximate to xo if and only if xe[a@, fol. 

Definition 2-2 Let R” be n-dimensional real vector space, U be an n-dimensional 
subspace of R”, xo=(X1,, X2- Xng) EU, and CcU be a "circular" region with center xo, 
which is called an approximation domain of xo. For Vx=(x1, x2,..., Xn)EU, say x is 
approximate to xo if and only if xEC. 

Definition 2-3 Let R” be n-dimensional real vector space, U=[a1, b;]<[d2, b2]x...x[dn, 
byJCR", and xo=(%X1,, X2)50-+. Xn, JEU. For Vx=(x1, X2,..., x)EU, say x is strictly 
approximate to xo if and only if the components x1, x,..., x, Of x are approximate to the 
components X1,, X255--+» Xno of xo, respectively, i.e., xelæ, Ailcla, bå (La, p is 
approximation domain of Xig> i=1,2, ... , n), and the "square" region [œ, (1]x[a@2, fo] x 
..-X [@n 8,]CU is called strict approximation domain of xo. 

In contrast to the strict approximation and strict approximation domain in Definition 
2-3, we call the approximation and approximation domain defined in Definition 2-2 to be 
ordinary approximation and ordinary approximation domain. It can be seen that the 
ordinary approximation is the approximation at the level of points (vectors), while the strict 
approximation is the approximation at the level of the coordinates of points (components of 
vectors). The relation between the ordinary 
the 
approximation domain of the same point xo 


approximation domain and strict 
is shown in Figure 2-1. The illustration also 


shows visually the relationship between the 


ordinary approximation and the strict w 
approximation. 
We know that the approximate Figure 2-1 Illustration of the relation between 


ordinary approximation domain and strict 


evaluation of a function is to obtain the ae : 
approximation domain 


approximate value of f( x’) in the condition ; : ; f 
Where the circular region is the ordinary 


of a pair (Xo, Yo) of corresponding values of 
function y=f(x) (x=(x1, X2,..., X,)) and x’ 
approximate to Xo being known. By the 


approximation domain of point Xo, and the 
square region is its strict approximation domain. 


axiom of approximation above, when x’ is approximate to xo, f(x’)=y is approximate to 
flxo)=yo. But approximation is a flexible concept (i.e., flexible linguistic value). In order to 
characterize precisely these two approximation relations and the relation between them, we 
introduce the measure of approximation-degree. 

Definition 2-4 Let R be real number field, xoe[a, b]CR, and [a, folcla, b] be the 
approximation domain of Xo. Set 


xQ-x x-a 


1—-——_=-_ .,, xe[a%, xo] 


x9— ao xoa 
Ay, (x)= (2-1) 


= x—p; 
1-a Tna “Eo A 


to be called the approximation-degree of x to xo. Where xo—Q=rl is called left approximation 
radius of xo, So—Xo=rr is called right approximation radius of xo. We call the function relation 
defined by the Formula (2-1) to be the approximation-degree function of xo. 

It can be seen that the range of approximation-degree function A,,(x) is [0, 1] and it is 


also reversible. In fact, it's easy to obtain that 


d (xo A)+Q 
Ax (x)= { (2-2) 


d (xo-po)+ Bo 
where d, €[0, 1] is the approximation-degree of x to xo. 
Definition 2-5 Let R” be n-dimensional real vector space, U be an n-dimensional 
subspace of R”, Xo=(X1,, X29»: Xn,)€U, and CCU be an approximation domain of xo. Set 


= Apacs xi)" 
Ag (=| 1 Qa, a). tn EC (2-3) 


to be called the approximation-degree of x to xo. Where r is the radius of C, and which is 
called the approximation radius of xo. 

We further consider the dependence relation between approximation-degree A, (y) and 
Ax, (x) under the constraint of function relation y= f(x). By Approximation Axiom, the closer 
x is to Xo, the closer y is to yo, thus the closer A, (y) is to Ax (x). Then when the 
corresponding approximation domains are small, we can simply agree that the two are equal, 
that is, set 


Ay, (W)=Ax, (x) 


This equation is equivalent to transmitting the approximation-degree of the independent 


variable x to xo to the corresponding function value y as the approximation-degree of y to yo. 
3 Approximate Evaluation of Functions Based on Approximation-Degree 


3.1 Approximate Evaluation of Univariate Functions Based on Approximation-Degree 


Let R be real number field, U=[a, b]<R, V=[c, d]cR, y=f(x) be a continuous and non-chaotic 
function relation from U to V. In the condition that a pair (xo, yo) of corresponding values of 
function y=f(x) and x’ approximate to x9 being known, find the approximate value of f(x’). 
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Let the approximation domain of xo be [a), bı]c[a, b]. According to the definition of the 
approximation-degree function above, the approximation-degree function of xo is 


=F, xe[a, Xo] 

Ax (x)= (3-1) 
x—bı 
EN, xe[xo, bı] 


And let the approximation domain of yo is [c1, d:]c[c, d], the approximation-degree function 


of yo is 


yey 
your Vale yo] 
Ay, = (3-2) 
y—dı 
ET yelyo, dı] 


The inverse function of A, (y) that is not difficult to get is 


dy(Yo — c1)+cı, dye[0, 1] 


Ay, (vy) = l (3-3) 
dyo — dı)+dı, d,€[0, 1] 


where d, is the approximation-degree of y to yo. 
The graphs of the approximation-degree functions A,,(x) and A, (y) are shown in 


Figure 3-1 and Figure 3-2, respectively. 


i Ax, (x) a Ay, 0) 
1 1 
0 0 
a Xo bi * a y d Y 
Figure 3-1 Graph of Figure 3-2 Graph of 
approximation-degree function Ay, (x) approximation-degree function Ay, (y) 


Now we find the approximation-degree A,, (x°), and then set Ay, = Ax, (x’) (that is, 
transmitting the approximation-degree of x’ to x9 to y’); Further, we derive the required 
approximate value of y' from approximation-degree A,,(y') and inverse function Ay, (yy 
of approximation-degree function of Ay, (y). 

It can be seen that the inverse function Ay, (y)' of Ay, (y) is a piecewise function, 
which has two parallel expressions. Thus, substituting approximation-degree Ay, (y’)=d into 
Ay, (y)', we can get two y-values (y; and y2). That is to say, from one x’ we obtain two 
y-values (as shown in Figure 3-3). Then, which y-value should be chosen as the desired 


approximate value of function ’? 


xX XxX x yı Yo X2 y 


Figure 3-3 One x’ is related to two y-values 


Obviously, the desired y-value is related to the orientation of x’ relative to xo and the 
trend (i.e., being increasing, decreasing, or a constant) of f(x) near xo. Thus, when the 
derivative f(x) is not known, we can consider whether there is a point x* on the x’ side near 
the point x», whose corresponding value of function, f(x*)= y*, is known. If there is such a 
point x*, we can estimate the trend of f(x) between the xo and x* by utilizing the size relation 
between the corresponding y* and yo, and then determine the choice of y-value. In fact, when 
x* <x’<Xo, if y*<yo, which then shows that the general trend of function f(x) is increasing on 
the sub interval (x*, xo), thus the y;, i.e., that y-value less than yo, should be chosen; while if 
y*>yo, which then shows that the general trend of function f(x) is decreasing on the 
subinterval (x*, xo), thus the ys, i.e., that y-value larger than yo, should be chosen. Similarly, 
when xo<x ’<x*, if yo<y*, which then shows that the general trend of function f(x) is increasing 
on the subinterval (xo, x*), thus the y should be chosen; while if yo> y*, which then shows 
that the general trend of function f(x) is decreasing on the subinterval (xo, x*), thus the yı 
should be chosen. But, if y*=yo, which then shows that the function f(x) is likely a constant on 
the subintervals (x*, xo) or (xo, x*), thus y=yọ can be taken no matter x*<x ’<xp or xp<x’<x*. 


3.2 Approximate [Evaluation of Multivariate Functions Based on 
Approximation-Degree 


Let's take the function of two variables as an example to discuss this problem. 

Let z=f(x, y) be a function (relation) from [a;, bı]x[az, b2] to [c, d]. Suppose a pair of 
corresponding values, ((xo, yo), Zo), Of function z=f(x, y) are known, and point (x’, y’) is 
approximate to point (xo, yo). In the case that expression of function f(x, y) is unknown or not 
used, find the approximate value of f(x’, y’). 

This problem is similar to the problem of function of one variable above, so we extend 
the above treatment to function of two variables, that is, we can find the approximate value of 
the function at point (x’, y’) from the approximation-degree of point (x’, y’) to point (xo, yo). 
Concretely, firstly get the approximation-degree Ac,,y,)(xy') , then set 
Az (Z)=A (xoy) (1y); and then, substitute it into inverse function A,, (z)! of A,, (z). Here, 
there will also be two z-values. Thus, we can consider whether there is an adjacent point (x*, 
y*) of point (xo, yo) in the direction determined by points (xo, yo) and (x’, y’), whose function 
value is known. If there exists such a point (x*, y*) (as shown in Figure 3-4), then we can 
determine the choice of z-value by taking f(x*,y*) as a reference. 

However, although the technique is feasible in theory, it has precondition, it requires that 
the points (x*,y*), (xo, yo) and (x’, vy’) must be exactly on a straight line. So, the method has 
some limitations for multivariate functions. Then, is there a better way? 
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By the definition of strict approximation 


(see Definition 2-3), point (x, yò is (xo, Yo) 
approximate to point (xo, yo) is equivalent to x’ 

is approximate to xo and y’ is approximate to yo. we 
Thus, we can find the approximate values z, (x**, y*) 


and z, of function f(x, yo) and f(Xo, y) at x’ and y’, 
respectively. Thus, we further imagine that if 
there is respectively an adjacent point of point Figure 3-4 Utilizing the value of the 
(xo, yo) in the x-direction and y-direction of the function at point («*,y*) to determine the 
point (xo, Yo), which is respectively (x*, yo) and approximate value of f(x’, y’) 

(xo, y*), as shown in Figure 3-5, and the 
corresponding function values f(x*, yo) and 
fo, y*) are known, then the approximate (x*, yo) Œ’, yo) a, 
value z, of fix’, yo) can be got by utilizing fix*, 
yo), and the approximate value z, of Axo, y’) O--| o y) 
can be got by utilizing Axo, y*), just like that 
of the previous univariate function. Thus, we 
firstly find separately the (Xo, »*) 


approximation-degree A% (x) and Ay, 0. 
then set Az (Z) =Ax, (x), and 4A, (z) = 
Ay, (y’); and then, substitute them separately 


Figure 3-5 Utilizing the values of the 
function at points (x*, yo) and (xo, y*) to 


. . . -1 
into inverse function Azo (z) of Azo (z) and determine the approximate values of f(x’, yo) 


get two pairs of candidate z-values, then . 4 Axo, Y’), respectively 

taking f(x*, yo) and f(x, y*) as reference, 

choose z, and z, from respective candidate values, respectively (as shown in Figure 3-5). 
Having got the approximate values z, and z,, how can we further obtain the approximate 

value z we need? Can z, and z, be synthesized into one z? For example, take the average value 

or weighted sum of z, and z, as z. But numerical experiments show that although the average 

value or weighted sum can be taken as the synthesis value, but the effect is generally not very 

ideal. However, fortunately, on the basis of the average value, by analyzing further, we have 

found a better way to synthesize z, and zy. 
Let z)=(z,+z,)/2 be the average 

of zx and z,. It can be seen from 

Figure 3-6 that z; can actually be 

viewed as an approximate value of 

f(x, y) at midpoint (denoted by (x, 

yi)) between (x’, yo) and (xo, y’). 

We can see from the figure that 


Z1<Z, Le., the varying trend of 

function values from to is a . 
; "A K Figure 3-6 Illustration 1 of synthesizing zx and zy into z 

decreasing. Set Zo-zi=c; (cı is the 

length of segment BC in Figure 3-7(a)), then z;=zo—c). Then, according to the varying trend of 

function values from Zp to zı (i.e., the slope of segment AB in Figure 3-7(a)), also, taking into 


account that point (xı, yı) is just the midpoint of segment joining points (x’, y’) and (xo, Yo), 
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that is, ||(x5 y) — (xo, yo)ll=20]| (x1, y1) — (xo, ¥o)I|, so we infer that the approximate value 
of function at point (x’, y’) can be zp—2c) (as shown in Figure 3-7(a)). Thus, it follows that 


Z=2o—-2€1=Zo—2(Zo—Z1)=Zo—-2[Zo-(Zx-Zy/2] = Zx+Zy —Zo 


Of course, zı may also be greater than zo or equal to Zo. If z,>zo, then the varying trend of 
function values from zo to zı is increasing (as shown in Figure 3-7(b)). Set zo—z1=c2 (c2 is the 
length of segment BC in Figure 3-7(b)), then z;=zo—c2. Then, according to the varying trend of 
function values from zo to zı (i.e., the slope of segment AB in Figure 3-7(b)), we infer that the 


value of function at point (x’, y’) can be zo—2c2. Thus, 


Z=Z0+2C2=Z0+2(Z1—20)=Z0+2[ (Zy+2Zy)/2—Zo] = Zet+Zy—Zo 


(a) (b) 


Figure 3-7 Illustration 2 of synthesizing zx and zy into z 


The third case: z;=Zo. This indicates that the values of function remain unchanged from zo 
to zı. Thus, we can take z=zo. And by zo=z1= (z,+2,)/2, it follows that 2z9=z,+z,. Thus, 


Z=229 —Z0= Zx+Zy —Zo 


In summary, we see that, no matter what relationship may be between the average of z, 
and z, and the zo, or no matter how the value of the function varies from zo to zı, the 
approximate value of the function at point (x’, y’) can always be taken as 


Z= Zx+Zy —Zo (3-4) 


Synthesizing the above analysis, we obtain a method for finding the approximate value 
of the function at the point (x’, y’) in the condition of that there are adjacent points (x*, yo) 
and (xo, y*) of point (xo, yo). That is, find firstly the approximate values of f(x’, yo) and fio, y `), 
Zx and z,, respectively, by f(xo, Yo)=Zo and approximation-degrees A, ED) and A, (y’) and 
by utilizing f(x*, yo) and f (xo, y*), and then synthesize zx and z, into a value z by Formula 
(3-4) as the approximate value of f(x’, y’). 

Now we see that here we separate actually the approximate evaluation of a function of 
two variables into the approximate evaluation of two functions of one variable, then, 
synthesize two obtained approximate values into a value as an approximate value of the 
original function of two variables. 

It is not difficult to see that the idea and method of “first separating then synthesizing” in 
the approximate evaluation of the functions of two variables above can also be extended to 
the approximate evaluation of the functions of three variables. That is to say, in the condition 


J 


of a pair ((%o, Yo, Zo), Uo) Of corresponding values of function of three variables, u= f(x, y, z), 
being known and point (x’, y’, z’) being approximate to point (Xo, yo, Zo), if there is 
respectively an adjacent point of the point (xo, yo, Zo) on the sides of corresponding 
coordinates x’, y’ and z’of point (x’, y’, z’) in the x-direction, y-direction and z-direction of the 
point (xo, Yo, Zo), which respectively is (x*, yo, Zo), (Xo, y*, Zo) and (xo, Yo, 2*) (as shown in 
Figure 3-8), and the corresponding values of function, f(x*, Yo, Zo), fo, y*, Zo) and f(%o, Yo, 2*), 
are known, then we can find firstly the approximate values of f(x’, yo, Zo), AXo, Y’, Zo) and f(xo, 
Yo, Z’), Ux, Uy and u, respectively, by f(xo, yo,Z0)=uo and approximation-degrees Axo (x", 
Ap (y) and A, (Z) and by utilizing f(x*,yo,z0), fxo, Y*,Z0) and fxo, yo, z*), and then 
synthesize the u,, u, and u, into a value u as the approximate value of f(x’, y’, z’). 


(xXo,Yo, 2*) 


(Xo0,Yo» 


! (x ’,Yo, Zo) 


(X0,¥0> Zo) 7 
(Xo, y : (x*,Yo, Zo) 


(Xo, Y*,Zo) 


Figure 3-8 Utilizing the values of the function at points 
(x*,y0,Z0),(X0,¥*,Z0) and (xo,Yo,Z*) to determine the approximate values 


of f(x’, Yo,zo), Axo, ¥’,Z0) and fxo, yo,z’), respectively 


From the synthesizing formula of approximate values of the function of two variables, 
Z=Zy4Zy-Zo=Z+Zy-Llezo, we infer that the synthesizing formula of approximate values of the 
function of three variables should be 


U= U,t+Uyt+U, —2Ug (3-5) 


In fact, when (x’, y’, z= (xo, Yo, Zo), Ux=Uy=u;=uo, While on the other hand, uo= 
UptUgt+Uo—2Uo. This indicates that for the especial point (xo, yo, Zo), the Formula (3-5) is correct. 
This indirectly verifies the rationality of the synthesizing formula (the three-dimensional AD 
interpolation example in Section 5.2 further confirms that the Formula (3-5) is reasonable). 

Based on the above analysis, generally, in the condition of a pair (Xo, yo) 
(Xo=(%1,.%2,>- re on a) of corresponding values of function of n variables, y=f(x) (x=(x1, X2, ... 5 
x,)) being known and point x ’=(x,’, x2’, ... , X,’) being approximate to point xo, if there is an 
adjacent point of the point x» on the side of corresponding coordinate of point x’ in the each 
coordinate direction of point xo, that is there exist points (x*, X29 >+ Ong ), (Xio ; 
Wi 52 T )avnny and (Mp X207: x,*), and the corresponding values of function are known, 
then we can find firstly the approximate values of f(x1’, X295- -Xngo) Xip: X25- Xng) > 5 
and fX10X20> Xn )s Vey» Yxz» ++» Yx,» respectively, by f(xo)=yo and approximation-degrees 
Ay, (%1'), Axg(*2') > -> Axo Xn") and by utilizing the values of function at these adjacent 
points, fri*, X%295-+-%Xng)> Jits Xn) JFX1o X20% Xn*), and then synthesize the 
Yxi» Yxz» Yx, into a value y as the approximate value of f(x’), and the formula of 


synthesizing these approximate values is 


n 


y=) y, -(n-1)yo (3-6) 
i=l 
In this way, we obtain an effective approach and method to find approximate values of 
multivariate functions. 
For convenience of narration, we may as well refer to the formula (3-6) as the 
Sum-Times-Difference formula. 


4 Interpolation Based on Approximation-Degree 


Let y=f(x) be a function (relation) from [a, b] to [c, d], and {(x;,y;)}7_, be a set of pairs of 
corresponding values of function y=f(x). Require to construct an interpolating function g(x) in 
the case that the expression of function f(x) is unknown or not used, such that g(x)= fx) G=1, 
2,..., n), and for other xe[a, b], g(x)* f(x). This is the usual interpolation problem. 

The readers may have noticed that the method finding functional approximate value 
based on approximation-degree above can also be used for interpolation. Below we discuss 
specifically this problem. 

Let a=x,<x2< , ... , < x,=b, then {x;}#_, is a set of interpolation base points (or nodes). 
We definite the approximation domain of x, is [x;, x2], the approximation domains of x; are 
[Xi-1, Xn] (i=2, 3, ... , n—1), and the approximation domain of x, is [x,-1, Xn], and definite the 
approximation-degree function of base point x; is 


x-X 
Ay, ŽE, xel, x] (4-1) 
17x2 
The approximation-degree functions of x; is 
X—Xj-1 
——,  xE[Xj_1,%; 
Xi—Xi—1 [ ii il 
Ay, (x)= (4-2) 
X—Xi+1 
——, xelx;, xi 
Xi—Xi+1 l i i] 


And the approximation-degree function of x, is 


Ay. @=——, xelxni Xn] (4-3) 


Xn—Xn-1 


The graphs of these functions are shown in Figure 4-1. 
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Figure 4-1 Graphs of the approximation-degree functions of numbers x1, X1,..., Xn 


Notes: 

e Here, the approximation domain of a number x; is determined not only by x; as before, 
but by two neighbors of x;, x;-; and x;,;.Therefore, two approximation radii rl= x; —x;, and rr= 
Xiz1-X; Of the number x; defined by Equation (4-2) are not necessarily equal. For the sake of 
distinction, we call the approximation domain determined by a number alone to be the 
absolute approximation domain of the number, and the approximation domain determined 
by its two neighbors to be the relative approximation domain of the number. In the absence 
of a special statement, the subsequent approximation domains refer to the relative 
approximation domains. 


Hat 


e When xe[x;_1, =), Ax; (x)<0.5 and Ay, _; (x)>0.5, this shows that x is closer to 


amin 


Xi-1; similarly, when xe( »Xi41], Ax; (x)<0.5 and Az, ,, (x)>0.5, this shows that x is 


Xi+1 


Xi itxj- 1 Xi man 
= —— A 


closer to x1; and only when xe[— ], Ax; (x)20.5. If we make an appointment: a 


point x is practically approximate to the base point x; if and only if its approximation-degree 


Xi a H Xi Sinu 


Ax; (x)20.5, then interval [—= ] is the practical domain of approximation-degree 


function Ax, (x), and that is also, the practical approximation domain of x; (while interval 

[Xi1, Xis1] is the conceptual domain of A, (x) and conceptual approximation domain of 

x;). And then, in the sense of practical approximation, we have the following conclusion. 
Proposition 4-1 A point x is practically approximate to the base point x; (i.e. 


A,,,(x)20.5), if and only if x lies in the practical approximation domain (ee a) of 


the x;. 

We then define the approximation-degree functions of y; (=1, 2, ... , n) in the same 
principle and way. 

First, we point out that yi, yo, ... , Yn corresponding to xı, X2 , ...., Xn in order do not 
necessarily satisfy yi<yo< , ... , <Yn. In fact, the size relationship between y; and yj-; in the 
sample data may be y;<yi-1, y=y;i-1, OF y;>yi-1; the size relationship between y; and yj,; in the 
sample data may be yi<yi+ı, Yi=Yi+1, OF Yi>Vi+1 . We also noticed that when x lies between x;_1 
and x;, the corresponding y lies certainly between y;_; and y;; and when x lies between x; 
and x;41, the corresponding y lies certainly between y; and yj,;. Thus, we only need to 
consider separately four (semi) approximation-degree functions of number y; in the four cases 
Of Vi-1<Vis ViSVi-1, ViSViv1, AN Vin 1 <Vi.- 

(1) Suppose yi-ı<y;; We define the first sub expression of the approximation-degree 
function of y; in this case to be 


y-Yi _ _Yi-1 i i 7 
Ay. OV" an 1 "Yi-Vi- ma Yi-Vi-v Pein yil SS 2 


Its graph is shown in Figure 4-2(a). 
(2) Suppose y;<yi-ı.- We define the second sub expression of the approximation-degree 
function of y; in this case to be 


i 1 i- 
AO aie nme Vee Yea] o> 
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Its graph is shown in Figure 4-2(b). 

(3) Suppose y;<y;,;. We define the second sub expression of the approximation-degree 
function of y; in this case to be 

Ay, (y)= ey, yeli Yil (4-6) 

Its graph is shown in Figure 4-2(c). 

(4) Suppose yi+ı<yi We define the first sub expression of the approximation-degree 
function of y; in this case to be 

1 Yi+1 


sd Iiia t y= ; : É 
AnS a pya” Vi-Viev yeliv yil G 7) 


Its graph is shown in Figure 4-2(d). 


(a) (b) 


Yi+1 Yi y 


(c) (d) 
Figure 4-2 Approximation-degree functions of number y; under different conditions 


Obviously, in the four expressions above, (4-4)=(4-5) and (4-6)=(4-7). That is to say, no 
matter y;1<y; or yi<yii, When y lies between yı and y, the expression of 


approximation-degree function of y; is always 


Y-Yi-1 1 Yi-1 
A.. SY S 
xO) Yi-Vi-1 VigVi-1 y Yi-Vi-1 


and no matter yjij<y; or yi<yiz1, when y lies between y; and yiu, the expression of 


approximation-degree function of y; is always 


Y—Yi+1 1 Yi+1 
A.. aA eg E 
xO) Yi-Vit1 YiTYi+1 Yi™Yi+1 


Thus, the above 4 functional expressions can be reduced as two expressions. 


Further, it can be seen that the coefficients of the first degree terms, i 
Yi-Vi-1 Vi-Vi+1 


in these two expressions are just the slopes of the corresponding straight lines. So they being 


positive or negative decide the corresponding function to be increasing or decreasing. Thus, 


itXj-1 
3 x;] and 


function y=f(x) being incremented or decremented on the subinterval == 
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Xi min 


[X;, ] are respectively determined by the signs of differences y;—y;_, and y;—y;41. 


This means that the numerical calculation takes the place of the logical judgment for the trend 
un i—1 XitXi+1 ] 


of function f(x) on the subintervals [ > 


A xi] and [xi 


And then, we get the inverse expressions of these two functional expressions, 
respectively 
Ay, (y) "=v ti (4-8) 
and 


Ay, (y V =d HV) +Y (4-9) 


Here d=Ay, (y)€[0, 1] is the approximation-degree of y to y;. 


Xi ma L Xi as 


Now, we transmit the approximation-degree of x (xe[— ) to x; to the 


corresponding y, that is, set 
dy=d,=A,, (x) 


et 


Also, considering that on the subinterval [ , Xi], the interpolated function y=f(x) may 


being increasing, decreasing, or a constant; while when y=f(x) is increasing, certainly y;1<y; 


ai 1 


so the corresponding ye[ , yil; when y=f(x) is decreasing, certainly y;<y;-ı, so the 


ae 


corresponding yely;, ~ ]; and when y=f(x) is a constant, y=y;1, so the corresponding 


yet nii 


1). Thus, when xe 


, yi] as well as ye[y;, , Xi], the corresponding 


y=yie[-——— 


X—Xj-41 


y; and y;, are adjacent, and here A, Om thus, the expression of the corresponding 


inverse function Ay, (vy) in Equation (4-8) becomes 


X-Xji-1 


-1_ 
Ay, O) = 0 iit 
_ WT Yi y p Fin Fi Vi 
Xi—Xi—1 Xi—Xi—1 
namely 
Yi—Yi—1 XiYi—-1—Xi—1Yi XitXi—1 
y y p SE yeH, x] (4-10) 

Xi—Xi—1 Xi—Xi—1 


Xi Hi 


Similarly, on the subinterval [x;,— ], the function y=f(x) may being increasing, 


decreasing, or a constant; while when y=f(x) is increasing, certainly y;<y,;, so the 


corresponding yel yj, zti y, when y=fx) is decreasing, certainly Yi+ı < y; so the 


yet 


corresponding ye[—~—, yi]; and when y=f(x) is a constant, y;=y;+ı, so the corresponding 


H nii 


Sey 


, Yi] as well as yeLy;, |. Thus, when xe[x;, , the corresponding y; 


y=yie(H 
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and yj; are adjacent, and here A, (x)= SS thus, the expression of the corresponding 
i+14%i 


inverse function Ay, (vy) in Equation (4-9) becomes 


-j Xj41-* 
Ay, O) =a q Om) +Y 
l l 
Yi Yir x+ XiVit1—*Xi+1 Vi 
Xi—Xi+1 Xi—Xi+1 
namely 
i—Yi XiYyi+17 Xi i xi+x; 
7i Vitd y 4 iYit+1 cy. xE[X;, i it (4-11) 
Xi—Xi+1 Xi—Xi+1 2 
; ; aT ; i+yi— 
Thus, with the two expressions, we can obtain directly the corresponding ye pe yi] 
YitYi-1 XitXi—1 a : Vit Vit 
or [y;, ne from xe[———., x;], and obtain directly the corresponding ye[y;, 7g 
it+Vi xitxi ‘ : 
or ee y] from xeE[x;, — ea but do not need to get approximation-degree and to 


evaluate inverse function. 

Actually, Equations (4-10) and (4-11) are two interpolating formulas. Thus, using these 
two interpolation formulas, we can solve the interpolation problem mentioned above. In fact, 
for the evaluated point xe[a, b], we only need to seek firstly the approximation domain the 
point x lies (This is equivalent to seeking nearest base point x; (i¢{1,2,...,n})), then to 
construct and use the corresponding interpolating formula according to the position of x 
relative to x; to do interpolation. 

In this way, we actually derive an interpolation method by the approximate evaluation of 
function based on  approximation-degree. We call this method to be the 
approximation-degree-based interpolation, or AD interpolation for short. 

Specifically, the practice of AD interpolation is: take base points (or, nodes) a= 
X1,X2,.+,X,=b as points of view, according to base points and their approximation domains 


xit? 
2 


to partition interval [a, b]=[x1, x,,] into 2n—2 subintervals as shown in Figure 4-3, [x4, 


|, 


X14 +xX2 
2 


Xat e E T F S . . . n 
2 5 tka x,,], as interpolation intervals; then, for evaluated point 


[ 


,X2), [X2, 


Xi +X j-1 


xe[ , Xi] do interpolating with Formula (4-10), for evaluated point xe[x;, ty do 


interpolating with Formula (4-11). Since each specific interpolating formula implies the trend 
of interpolated function y=f(x) on the corresponding subintervals, therefore, there is no much 
error between the obtained approximate value and the expected value, and there would not 
occur the case that two y-values are got from an x. 
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ag x1 +x2 7 X2 +x3 4 x3 +x4 Xn—2+tXn—1 Xn—1 Xn 
ua 2 2 2 oe i 2 n-1 2 
Figure 4-3 Illustration of interpolation intervals partitioned by base points and their 


approximation domains in one-dimensional interpolation 


Where the solid lines indicate the boundaries between approximation domains of base 


points x;, and a dotted line further divides an approximation domain into two parts 


Example 4.1 Use AD interpolation to do interpolating to functions y=sin x and y=cos x, 

respectively. 

We take sampled-data points of y=sin x, (x; y;), as follows: 

x; Oxa, 0.1 <a, 0.22, 0.3xz, ... , 1.92, 2x7. 

yg sin(x). 

and take evaluated points x: 0x, 0.02xz, 0.04xz, 0.06xa, ... , 1.98xz, 2xT. 
By using our AD interpolation to do interpolating (programmed with MATLAB), the 
corresponding y-values obtained are as follows: 

0 0.0624 0.1249 0.1873 0.2497 0.3118 0.368 1 0.4245 0.4808 
0.5371 0.5923 0.6369 0.6816 0.7263 0.7710 0.8133 0.8420 0.8707 
0.8994 0.9281 0.9530 0.9629 0.9728 0.9827 0.9926 0.9975 0.9876 
0.9778 0.9679 0.9580 0.9424 0.9138 0.8851 0.8564 0.8277 0.7934 
0.7487 0.7040 0.6593 0.6146 0.5653 0.5089 0.4526 0.3963 0.3400 
0.2809 0.2185 0.1561 0.0936 0.0312 —0.0312 -0.0936 -0.1561 —0.2185 


—0.2809 —-0.3400 —0.3963 —0.4526 —0.5089 —0.5653 —0.6146 —0.6593 
—0.7040 -0.7487 —0.7934 —-0.8277 —0.8564 -—-0.8851 -0.9138 -0.9424 —-0.9580 
—0.9679 —0.9778 —0.9876 —0.9975 —0.9926 —0.9827 —0.9728 —0.9629 


—0.9530 -0.9281 -0.8994 -0.8707 -0.8420  -0.8133 -0.7710 -0.7263 -0.6816 
—0.6369 -0.5923 -0.5371 -0.4808 -0.4245 -0.3681 -0.3118 -0.2497 -0.1873 
-0.1249 -0.0624 0.0000 
The interpolation effect is shown in Figure 4-4(a). 

To function y=cos x, we do interpolating with the same method and data, the 
interpolation effect is shown in Figure 4-4(b). 
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(a) (b) 


Figure 4-4 The effect drawings of AD interpolation to functions y=sin x and y=cos x 


Where the blue circles, “o”, in the graph indicate the sample data points, and the red 


crosses, "+", indicate the points obtained by AD interpolation. 
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Accuracy analysis: 
sin x: maximum error 0.0121, minimum error 0, average error 0.0052 
cos x: maximum error 0.0122, minimum error 0, average error 0.0052 
Actually, it can be seen from the interpolation formulas that one-dimensional AD 
interpolation is a local linear interpolation. That is to say that this one-dimensional AD 
interpolation has the effect of “reaching the same goal from different routes” with usual local 
linear interpolation. 


5 Multidimensional and High-Dimensional Interpolations Based on 
Approximation-Degree 


5.1 Two-Dimensional AD Interpolation 


Suppose some data points of function of two variables, z=f(x, y), (x, y), Zj) G=1,2,....n, 
J=1,2,...,m), are known, and the base points (x; y;) are regularly distributed, that is, they 
satisfy a=x\<x2< , ... , SXn=by and ay = yi<yo< , ... , <Vm=b,, and can also form an mxn 
matrix: 

(1, Yi), as Y1), +++ » ns Y1) 

(X1, y2), (2, Y2), +++ » Xn, Y2) 


(xı, Ym), a, Ym), e... (Xns Ym) 


Correspondingly, they form a Cartesian grid on the x-y plane (as shown in Figure 5-1). We 
consider corresponding interpolation problem. 


X4+X2 x X2+X3 m X3+x4 
2 2 2 a 02 4 


X41 


Figure 5-1 Illustration of interpolation regions partitioned by base points and their 


approximation domains in two-dimensional AD interpolation 


Where the regions enclosed by solid lines are the approximation domains of base points 


(x;, y;), and the dotted lines further divide an approximation domain into two or four parts 


For this kind of two-dimensional interpolation (actually, this kind of interpolation is also 
usual two-dimensional interpolation), we use firstly the idea and technique similar to AD 
interpolation in Section 4 above, that is, take base points (x1, y1), ... 5 (Xn Ym) as points of view, 
according to base points and their approximation domains to partition region U=[a,, b,|x[ay,, 


b,]J=[ x1, Xn ]X[ Y1, Ym ] into (2n—2)x(2m—2) sub-regions as shown in Figure 5-1 as 
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interpolation regions; then, for evaluated point (x, y)€U, seek the sub-region it lies (This is 
equivalent to seek the corresponding nearest base point (x, y,) Ge {1,2,....2}, je{1,2,...,.m}) 
and determine the position of (x, y) relative to (x; y,)), and then to do interpolation. 

However, for the interpolation on these two-dimensional sub-regions, we have not 
derived the corresponding interpolating formulas. Since these base points are regularly 
distributed, that is, each base point has the adjacent base point(s) in the direction of its 
coordinates, so we use the technique of “first separating then synthesizing” described in 
Section 3.2 above to get the interpolated value at point (x, y). Specifically, first separate 
interpolating computation of evaluated point (x, y) with corresponding nearest base point (x; 
y;) into interpolating computation of x with x; and interpolating computation of y with y,, and 
getting respectively the corresponding approximate values z, and z,; then synthesize z, and z, 
into a value z as the interpolated value at point (x, y). 

Obviously, the two interpolating computations obtained from the separation are already 
one-dimensional interpolating computations, so previous two interpolating formulas can be 
used. But, here the corresponding interpolating formulas are with respect to x and y, 
respectively, and which have the total of four formulas separated as two pairs: 


Zij —Zi—1 j XiZi—1,j —Xi-12ij XitXi—1 
ay p xe QE, ri 5-1 
{ xiki- Xi-Xj-1 [ 2. il (5-1) 
Zij —Zi+1,j XiZi+1,j Xi+1Zij xitxi 
= i t+ Diy t4i+1, i+ L xelx;, i +h) (5-2) 
Xi—Xi+1 Xi—Xi+1 2 
Zij —Zij—1 YjZij-17Yj-1Zij Yj +tYj—1 
a A el j] (5-3) 
| Vj-Vj-1 y Yj-yj-1 °? 7 2 Yj 
Zij —Zij+1 YjZij+17Yj+1Zij Yj+tyYj+1 
wey HAA, yey, SS] (5-4) 
YjYj+1 y ViVi” yew 2 


where zi;=fxi, y;). 

Of course, for an evaluated point (x, y)eU, we need only choose a formula with respect 
to x and a formula with respect to y. As to which two formulas will be used, it should be 
determined according to the subintervals that x and y lie respectively. For example, if 


e[ 


Xi tXxi-1 
2 


, x] and ye[y,, ama, then Formulas (5-1) and (5-4) should be respectively 
used. 

After getting corresponding z, and z, by two formulas used, it is must to synthesize them 
into a value z. From Formula (3-7) in Section 3.2, here the formula that synthesizes z, and z, 
should be 


Z= Zx+Zy —Zij (5-5) 


In this way, we have extended the AD interpolation to two-dimensional interpolation. 
From the Formula (5-5) it can be seen that this kind of two-dimensional AD 
interpolation of “first separating then synthesizing” is different from the usual bilinear 
interpolation, but its effect is comparable to that of the latter. 
Example 5.1 Use AD interpolation to do interpolating to function z=—-x°-y°. 
We take sampled-data points of z=-x y", (Œi; yj), Zi), as follows: 
xj: —20, -18, -16, ... ,—2, 0,2, ... , 16, 18, 20. 
y; —20,-18,-16, ... , —2, 0, 2, ... , 16, 18, 20. 
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zji ary 
and take evaluated points (x, y) as follows: 

x: —20, —20, 20, -19.5, —17.8, -18,-15.3, —12, —10.2, —10, —10, 0, 0, 10, 10, 5.6, 4.7, —3.4, -1.8, 
—2.3, -3.6, 1.2, -5.4, -15.6, —20, —20, —20, -18.3, 18.4, 17.5, 16.2, 14.5, 11.1, -5.4, -12.1, -8.5, -13.9, 
7.5, -7.8, —9.8, —12.4, -13.5, -14.6, -17.5, -17.8. 

y: —20, 20, —20, —19.5, -17.8,—5, —15.5, 2.5, -10.2, 10, —20, 0, —20, —20, —10, —15.3, —3.8, 
13.4, —2.8, —1.9, —5.6, —10.2, —6.5, 5.6, 0, —10, 10, 10.4, -18.1, —16.3, -14.4, -12.3, -6.3, -15.8, —8.2, 
—15.6, 0.9, 1.6, 3.2, 4.6, 6.6, 2.8, —0.9, 18.6, 13.2. 

Doing interpolation by using AD interpolation (programming with MATLAB), the 


corresponding z-values obtained are as follows: 
-800.0000 -800.0000 -800.0000 -762.0000 -634.4000 -350.0000 -476.0000 
—151.0000  -208.8000 -200.0000 -500.0000 0 -400.0000 -500.0000 
—200.0000 -267.0000 -37.8000 -192.8000 -12.4000 -9.6000 —45.6000 
-106.8000 -73.0000 -276.0000 -400.0000 -500.0000 -500.0000 -444.2000 
667.0000 -573.2000 -470.8000 -362.8000  -164.4000 -280.0000 -214.2000 
-317.0000 -195.2000 -60.2000  -72.4000 -1184000 -198.8000 -—191.8000 
-215.8000 -653.8000 -492.4000 

The interpolation effect is shown in Figure 5-2. 


Figure 5-2 The effect drawing of scattered AD interpolation to function z=—x"—y” 


Where the grid curve is the graph of function z=—x-y’, and the red circles indicate 


the points obtained by AD interpolation. 


Accuracy analysis: 
maximum error 1.8300, minimum error 0, average error 0.9222 


Example 5.2 Use AD interpolation to do interpolating to function z-i x? = ay. 


We take sampled-data points, ((x;, yj), zi), as follows: 
xi —5,—4, —3, -2, -1, 0, 1, 2, 3, 4, 5. 
yj. —5, -4, -3, -2, -1, 0, 1, 2, 3, 4, 5. 
1 1 
Zij: ee z AE 
and take evaluated points (x, y) as follows: 
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x=-5, —4.9, -4.8, -4.7, -4.6, -4.5, ... , 4.5, 4.6, 4.7, 4.8, 4.9, 5. 
y=x. 
Doing interpolation by using AD interpolation (programming with MATLAB, the 


corresponding z-values obtained are omitted), the interpolation effect is shown in Figure 5-3. 


(a) (b) 


Figure 5-3 The effect drawing of AD interpolation to function aon? = - y? 


Where the (a) is the functional graph before interpolating, and the (b) is the 
functional graph after interpolating 


Accuracy analysis: 
maximum error 0.0625, minimum error 0, average error 0.0216 
Example 5.3 Figure 5-4 below shows also an effect drawing of AD interpolation, the 
sampled-data (are omitted) are taken from the classic peaks function in MATLAB. 
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Figure 5-4 An illustration of the effect of AD interpolation 


Where (a) is the functional graph formed by sample data, and (b) is the functional graph 
obtained by AD interpolation. 


Accuracy analysis: 
maximum error 0.3177, minimum error 0, average error 0.0224 
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5.2 Three-Dimensional AD Interpolation 


Suppose some data points of function of three variables, u=f(x, y, z), (Œ; Yp Ze), Ui), 
(i={1,2,....n}, j={1,2,....m}, k={1,2,...,/}) are known, and the base points (x; y; Zx) are 
regularly distributed, that is, they satisfy a,=x)<x2< , ... , <m=D, , dy = Vi<yo< , «5 SYm=by 
and a, =Z|<Z2<,... , <z=b,, and can also form a mxnxl “three-dimensional matrix” (as shown 
in Figure 5-5), correspondingly, they form a Cartesian grid in the x-y-z space (as shown in 
Figure 5-6). We consider the corresponding interpolation problem. 


Bi i i (x1, yis zi), (X2, Yip zi), ... (Xn yi» 2° 
I P 


l 
l 
i 
/ d I 

/ I , 
e (X15 Yz Z), (X2, Yo Zi), «++ 5 Xnr V2, Z) i 
’ Pe I 
Bp eee eee age 
vil Ld l 
oy g I 
a“ i (Xi, Yi, 22), X2, Y1, Z2), «e+ » (Xn, Y1, 32) en Ym z) ! 
An neem nnn nee í i 
l l 


(X1, Yo Z): n Yis Zi) e > Ons Yis Z) 


(X1, Yz Z1), X2 Yz Z1), «++ An Y2 Z1) 


(x1, Yms zı), (X2, Yms zı), a (Xn; Yms zı) o 


Figure5-5 A “three-dimensional matrix” formed by known base points, (x;, Yj, Z) 


For the three-dimensional interpolation, we use still the technique similar to the 
two-dimensional interpolation above. Specifically speaking, it is that take base points (x, yı, 
Zi) «+» » (Xn Ym Zi) as points of view, according to these base points and their approximation 
domains to partition region U=[a,, b,]x[a,, b,]x[a,, b:] into (2n—-2)x (2m—2)x (2I-2) sub 
regions as shown in Figure 5-6 as interpolation regions; then separate the interpolating 
computation of evaluated point (x, y, z) with corresponding nearest base point (x; yj, Z) 
(ie {1,2,....n}, je {1,2,...,.m}, ke {1,2,...,/}) into the interpolating computation of x with x;, the 
interpolating computation of y with y,, and the interpolating computation of z with zg, and get 
the corresponding approximate values ux, uy, and u,, respectively; finally, synthesize u,, uty, and 
u, into a value u as the interpolated value at point (x, y, z). 
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X4+X2 x X24X3 x xX34+X4 x x 


X1 7 2 7 3 2 4 
Figure 5-6 An illustration of interpolation regions partitioned by base points and their 


approximation domains in three-dimensional AD interpolation 
Where the regions enclosed by solid lines are the approximation domains of base points (x;, yj, Zi), 


and the dotted lines further divide an approximation domain into two, four, or eight parts 


There were 4 formulas separated as 2 pairs in the above two-dimensional interpolation, 


for here three-dimensional interpolation we have 6 formulas separated as 3 pairs: 


Uijk —Ui-1 jk XjUj-1 jk —Xi-1Uijk Xj+Xj-1 
y= OE y xe ES, y] (5-6) 
Xj—Xj-1 Xj-Xi-1 2 
Uijk —Ui+1 jk XjUit1 jk —Xi4+1Uijk XitXi+1 
y= y p EE elx, 1H (5-7) 
Xi—Xi+1 Xi—Xi+1 2 
_Yijk Ui j—1,k YjUij-1,.k—Vj-1Uijk Yj+Yj-1 
u= y+ > yel 2 , yj] (5-8) 
Vji-Vj- Vj-Vj-1 
Uijk “Ui j+1,k VjUijtik—Vj+1Uijk Yj tYj+1 
u= y+ > yely;, ] (5-9) 
Yj—Yj+ Yj—Yj+1 2 
Uijk —Uij,k—-1 ZķUijk—-17Zk—1ťijk Zk tZķk-1 
=s y p H, zeA, 2x] (5-10) 
l Zk -Zk— Zk -Zk-1 2 
Wijk —Uij,k+1 ZU k+1—Zk+1Uijk ZktZk+1 
u= z + oe cr, E ] (5-11) 
Zk—Zk+ Zk—Zk+1 2 


where u= fxi, Yj, Z). 

Similar to two-dimensional AD interpolation, here for an evaluated point (x, y, z) €U, we 
need only choose a formula with respect to x, a formula with respect to y, and a formula with 
respect to z. As to which 3 formulas will be used, it should be determined by the subintervals 


—1 


; ; , +x; ajo 
that x, y and z lie respectively. For example, if xen, xil, yely;, wa, and 


zeka, Zęl], then Formulas (5-6), (5-9), and (5-10) should be respectively used. By 


Formula (3-12) in Section 3.2, here the formula that synthesizes u,, uy, and u, is 
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U= Uy+Uy+Uz —2Ujx (5-12) 
Example 5.4 Use AD interpolation to do interpolating to function of three variables, u= 


-x-y z. 
We take sampled-data points, ((x;, Yj, Zk), Wijk ), as follows: 
x; —3, —2.5, —2.0, ... , —0.5, 0, 0.5, ... , 2.0, 2.5, 3. 


y; =—4, -3.5, —3.0,..., —0.5, 0, 0.5,..., 3.0 , 3.5, 4. 
Zi —4, -3.5, -3.0,..., -0.5, 0, 0.5,..., 3.0, 3.5, 4. 


. = 2 2 2 
Wijk: U= Xj -Yj Zk 


and take evaluated points, (x, y, z) as follows: 


x: —3, -2.75, -2.50, ... , -0.75, 0, 0.25, ... , 2.50, 2.75, 3. 
y: —4, -3.80, -3.60, ... , -0.80, 0, 0.20, ... , 3.60 , 3.80, 4. 
z: —4, -3. 75, —3.50, ... , -0.75, 0, 0. 25, ... , 3.50, 3.75, 4. 


Doing AD interpolation (the corresponding u-values obtained are omitted), the interpolation 
effect (slice chart) is shown in Figure 5-7. 


Figure 5-7 The effect drawing (slice chart) of AD interpolation to function of three 


variables, u= -x-y =z’ 
Where the graph on the left is the graph before interpolating, and the graph on the right is 
the graph after interpolating 
Accuracy analysis: 
maximum error 0.1850, minimum error 0, average error 0.0993 
5.3 N-Dimensional AD Interpolation 


Generalizing the above two-dimensional and three-dimensional AD interpolations, we obtain 
a general n-dimensional AD interpolation method. 
Suppose some data points of function of n variables, y= f(x, X23, ... 5 Xn) 


(1X2) es Xn b Vij wk) (=1,2,...,r; JH1,2,...,5 ; k=1,2,...,t) are known, and the base points 
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(x1, X2; Xn) are regularly distributed, then the procedure of n-dimensional AD 


interpolation is: 
(1) Partition the n-dimensional interpolation space U=[a,, b,]x[a2, b2)x...x[a,, 
Dal=[ X1,,%*1,1x[%2,,X2,]x---X[%n,,Xn,] according to the approximation domains of all base 


points Xij K=p XZ, Xn) (=1,2,...,7; J=1,2,....8 ; K=1,2,...,6)3 


(2) For evaluated point x=(x1, X2, ... , X,)€U, look for the approximation domain it lies 
(Which is equivalent to seek its nearest base point x;j..« (ie {1,2,....7}, fe{1,2,....5}, 
ke {1,2,...,t})); 

(3) Separate the n-dimensional interpolating computation of point x with base point 


ar 


Xij, into n one-dimensional interpolating computations of x, with X41 ,> X2 with X2 js vers and 


Xn with Xn,, and do interpolating computations with the corresponding one-dimensional 
interpolation formulas, respectively; 

(4) Synthesize the n got approximate values y,,, Yy,,-+,Yx, into one value y as the 
interpolated value at point x. 


N-dimensional AD interpolation has n pairs of interpolating formulas as follows: 


Vij ..k ~Yi-1,,,k X1 Vi-1j,,k7*1,;_1V ij uk Xq,+%4,_ 
: a. ao a ej, (5-13) 
X4,—-X1,_ X4,—-%1,_ 2 
Vij wk “Viti j,k X1,Vit1j kX lj yi k X1, 4% 1344 
RPL Plt yy p SP TE ely, Ht] (5-14) 
li Alit i “lit 
Vij kY ij-1,„k X2Yij-1,nk7X2;_4Vij ..k X%2,+%2;_4 
ye Metal p p St eft xy] (5-15) 
2; *2;_ X2;—%2;_ 
Vij uk Vij +1,,k X22 Vijt1k7%2;44V ij uk X2;,4%2544 
aM Pat py p ET ela, Ht] (5-16) 
2i 7X24 X2 X24 
Yij k Yij,k-1 XniYijnk-17¥ni_1Yüj ...k Xn; +%Xnj_4 
ye MU yf ky eft xy] (5-17) 
ni~ 4ni ni tni 
Yij.kTYijnk+1 Xn Vij,,kt+1—%n it.) i „k Xn; +Xn i+1 
yer i yf S E, Xn E] (5-18) 
ni Anita ni Aniya 
where Yi _K=f(%1,; X2 vee Xn,)- 
And the synthesizing formula of the final interpolated value is 
n 
Y= Dy, -Dya (5-19) 


i=l 


Here yy, are the approximate values that got by one-dimensional AD interpolation 
computations of x; (i=1,2,...,n). 

The Formula (5-19) is the Sum-Times-Difference formulas in AD interpolation. 

It can be seen that in theory the dimensionality n of n-dimensional AD interpolation has 
no an upper limit. In this way, using the technique of “first separating then synthesizing” and 
the Sum-Times-Difference formula, we can also realize the high-dimensional or even 


super-high-dimensional AD interpolations with more than 3 dimensions. In fact, a program 
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for AD interpolation with more than 10 dimensions has been completed by the author, and the 
dimension can be easily increased (due to the space limit, the relevant content will be 
introduced in another article). 

Actually, if we substitute the interpolation formulas into the Sum-Times-Difference 
formula, it is not difficult to see that the Sum-Times-Difference formula is actually a linear 
combination of all coordinate components of point (x1, x2, ... , Xn). For example, we substitute 
two-dimensional interpolating formulas ge and ge y + 


Xj-Xj-1 Xi—Xi—1 Vj-Vj-1 


;Zij-1-Vj-12Zij . : : ; . 
HEATA 3 ae 1“! into corresponding Sum-Times-Difference formula z= z,+zy —z,, then it follows 
jv o 

that 


Zij —Zi—1,j Zij —Zij—1 XiZi—1j—Xi—-1Zij , Vj4ij-17Vj-12ij 
=U j ypt j + ( j j y iij j j i) 
Xj-Xj-1 Vj-Yj-1 Xj —Xj-1 Vj-Yj-1 


Therefore, Sum-Times-Difference formula can also be said to be the interpolation formula of 
multidimensional AD interpolation. 

Although AD interpolation is also a local linear interpolation, it is derived on the basis of 
approximate evaluation of functions based on the approximation-degree, so the base points 
involved in interpolation are related to the position of the evaluated point xeU relative to its 
nearest base point x;;__,. It can be seen that n-dimensional AD interpolation involves a total 
of 1+n base points, and when dimensionality is added by 1, the number of the base point 
involved will be added by only 1 (see Figures 3-7, 3-8 and 5-1, 5-5). Because the point x is 
only approximate to the base point x;;,..4, so the n approximate values, Yj; Yx» Yxn» 
obtained from one-dimensional interpolating computations got from separating are most 
affected by base point x;;,__;, and final interpolated value y is also most affected by x;;,...x- In 
other words, AD interpolation is dominated by the nearest base point x;;,__; of evaluated point 
x. 


6 Accuracy and Efficiency about AD Interpolation 


6.1 The Accuracy 


It can be seen from the one-dimensional interpolation formulas and the 
Sum-Times-Difference formula of multi-dimensional interpolation that AD interpolation is 
essentially a piecewise (block) linear interpolation. Therefore, the accuracy of AD 
interpolation should be on the level of linear interpolation. In fact, the data experiments show 
that: the accuracy of 1, 2, 3-D AD interpolations are about the same as the traditional linear 
interpolation, especially for the quadratic functions (such as previous examples 1, 2, 3), the 
accuracy of the two is exactly the same. In addition, we find that to the interpolated function 
y=f(x1, X2, ... , Xn) which is monotonous on each interpolation subinterval, in theory, the 
maximum error of AD interpolation is 
X1; = X2, +X2, Xn; +HXn a *1; i X2, +X2, „ni +Xy 
E=| foc iti i - iti, a js fH iti i 5 itt ae) 
Here, f(x) and f z are respectively the interpolated value and desired value at point 
X=(X1, X2, ... 5 Xp). 
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Actually, from the principle of AD interpolation, it is not difficult to see that the accuracy 
of AD interpolation is related to two factors: one is the size of the approximation domain of 
correlative base point and the other is the approximation-degree of the evaluated point x to the 
corresponding nearest base point. Obviously, the narrower the approximation domain 
(including the relative approximation domain and the absolute approximation domain) is, the 
higher the interpolation accuracy is; And when an approximation domain has be determined, 
the higher the approximation-degree of the point x to the corresponding base point is, the 
higher the accuracy of the corresponding interpolated value is. 

Thus, for the regularly distributed base points, as long as there are enough data points, 
that is, enough samples, the required any accuracy can be achieved by using AD interpolation. 


6.2 The Efficiency 


Viewed from the principle, the difference between AD interpolation and traditional linear 
interpolation is mainly in multidimensional interpolations. In fact, in multidimensional AD 
interpolation, by using a technique of "first separating and then synthesizing", an 
n-dimensional AD interpolation computation is separated into n one-dimensional AD 
interpolation computations to be done respectively; then the got results are synthesized into a 
value by Sum-Times-Difference formula as the result value of the n-dimensional interpolation. 
In this way, one-time n-dimensional AD interpolation computation becomes into n-times 
one-dimensional AD interpolation computations and one-time synthesizing computation. 
From the Sum-Times-Difference formula, all of multidimensional AD interpolations have 
only two layers of computation: one-dimensional AD interpolation computations and 
synthesizing computation (as shown in Figure 6-1 (a)), and when dimensionality is added by 
1, only one-time one-dimensional AD interpolation computation and one-time addition 
operation are added in the corresponding AD interpolation (as shown in Figure 6-1 (b)). That 
is to say, the two-layer computing structure of AD interpolation is always unchanged, but the 
number of terms of computation in lower layer will increase with the increase of data 
dimension. In this way, the times of the formula (including interpolation formulas and 
Sum-Times-Difference formula) computation in n(n>1)-dimensional AD interpolation is n+1. 
Starting from the two dimensions, the times of formula computation in the AD interpolation 
are in order: 3, 4, 5, ..., i.e., a sequence of equal difference numbers with a difference of 1, 
and the general term is: n+l. We denote the times of the formula computation in 
n(n>1)-dimensional AD interpolation as N,,, and the recurrence formula is: N,.,;= N,+1. From 
this we see that the computation load of multidimensional AD interpolation is much smaller 
than that of traditional multidimensional linear interpolation. 


(a) (b) 
Figure 6-1 The hierarchical structure of formula computation in 


multidimensional AD interpolation 


Where I,, denote the one-dimensional AD interpolating computations of 


variable x; (i=1.2.....2.n+1) 
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In addition, since n one-dimensional AD interpolation computations got from separation 
are independent of each other, so the parallel computation can be used. Thus, the efficiency of 
n-dimensional AD interpolation is almost equal to that of one-dimensional AD interpolation 

From the above mentioned, we can see that the high efficiency of AD interpolation is 
obvious, and the higher the dimension is, the higher the efficiency is. 


7 Summary and Prospect 


Inspired by the approximate evaluation method of flexible linguistic functions in reference [1], 
in this paper we introduce the measure of approximate-degree and the concept of 
approximate-degree function between numerical values, and according to "approximation 
axiom" and by computing approximation-degree, transferring approximation-degree, 
evaluating inverse function, and choosing approximate value, we realize the approximate 
evaluation of a numerical function; Especially, for multivariate functions, in the sense of 
"strict approximation", we use the technique of “first separating then synthesizing” and find 
and present a formula of "Sum-Times-Difference", and achieve corresponding approximate 
evaluation more skillfully. On this basis, we further derive a set of interpolation formulas for 
regularly distributed base points, thus an interpolation method based on approximation-degree, 
approximation-degree-based interpolation, i.e., AD interpolation, has be developed; and then, 
we extend one-dimensional AD interpolation to high-dimensional interpolation. As a test of 
AD interpolation method, we apply AD interpolation to the actual functions and obtain 
satisfactory results: the one-dimensional AD interpolation has the effect of “reaching the same 
goal from different routes” with usual piecewise linear interpolation, the two-dimensional and 
three-dimensional AD interpolations are comparable to the usual bilinear and trilinear 
interpolations; and then, the high-dimensional AD interpolations with more than 3 dimensions 
also have passed the data test. 

Actually, AD interpolation is also a kind of linear interpolation in essence, its accuracy 
is about the same as that of traditional (piecewise) linear interpolation, but its principle is 
different from that of the latter. One-dimensional AD interpolation is done directly by using 
corresponding interpolation formulas; n(m>1)-dimensional AD interpolation is firstly 
separated into n one-dimensional AD interpolation computations which can be processed in 
parallel to do respectively, and then the results got are synthesized by Sum-Times-Difference 
formula into a value as the result value of the n-dimensional interpolation. AD interpolation is 
dominated by a nearest base point and assisted by adjacent base point(s) of the nearest base 
point, n-dimensional AD interpolation computation involves 1+n base points, which means 
that when dimensionality is added by 1, only one correlative base point will be added. The 
n-dimensional AD interpolation has two layers of computation: first, parallel one-dimensional 
AD interpolation computations, then, synthesizing computation, and when dimensionality is 
added by 1, only one term of one-dimensional AD interpolation computation and one-time 
addition operation are added. In this way, compared with the traditional linear interpolation, 
the computation load of the AD interpolation is much smaller, if the parallel processing is also 
used, the efficiency of n-dimensional AD interpolation is almost the same as that of the 
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one-dimensional AD interpolation. In addition, the dimensionality of the multi-dimensional 
AD interpolation is theoretically unlimited (limited only by the associated computational 
resources). 

In a word, AD interpolation has the advantages of novel idea, unique method, simple 
computation, adjustable accuracy, high efficiency, and unlimited dimension. This opens up a 
new way for data fitting and function approximation, and especially, it starts a feasible and 
convenient approach and provides an effective method for high-dimensional and 
super-high-dimensional interpolations. 

This paper expounds the basic principle of AD interpolation. On the basis, we can 
further carry out the following work: 

e Applying AD interpolation to practical high-dimensional and super-high dimensional 
interpolation problems 

e AD interpolation with irregular distributed base points 

e AD interpolation for vector-valued functions 

e AD interpolation with derivative (partial derivative and directional derivative) 

e Using AD interpolation to realize a function approximator 

It can be seen that the approximator using AD interpolation is always interpretable, and 


won’t occur “the dilemma between precision and interpretability” ™ 


encountered by the 
function approximator based on fuzzy technology. 
e Introducing AD interpolation into machine learning to realize a new instance-based 


learning method 


AD interpolation learning 
Actually, for these further topics, the author has done some research and obtained some 
preliminary results. Due to space constraints, we will discuss them separately. 
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